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ABSTRACT 

We have used the VLA to detect emission from the supergiant VY CMa at ra- 
dio wavelengths and have constructed 3000-4500 K isothermal outer atmospheres 
constrained by the data. These models produce a radio photosphere at 1.5-2 R*. 
An extrapolation of the model can account for the observed total mass-loss rate 
of the star. We also present mid-infrared imaging of the supergiant that suggests 
warm dust is extended in the same direction as the near-infrared reflection neb- 
ula around VY CMa. The origin of the asymmetries in the outflow remains an 
unsolved problem. 

Subject headings: infrared: stars — radio continuum: stars — stars: atmospheres 
- stars: individual (VY Canis Majoris) — stars: mass loss — supergiants 



1. Introduction 

Massive stars (M iniUa i > 8-10 M ) conclude their stellar evolution as red supergiants 
before exploding as Type II supernovae. VY CMa (HD 58061; M5ela) is one of the best 
studied supergiants. At a distance of 1.5 kpc (Lada & Reid 1978), VY CMa maintains 
a luminosity of L* = 2-5 x 10 5 L (Jura & Kleinmann 1990; Le Sidaner & Le Bertre 
1996; Monnier et al. 1999), suggesting it is a massive star (M > 15 M ) that will explode 
as a supernovae within 10 5 yr (Hirschi et al. 2004). It is losing mass at a rate ~ 2-4 x 
10~ 4 M yr -1 (Danchi et al. 1994; Monnier et al. 1999) which produces the optically thick, 
dusty envelope enshrouding the star, thus hiding many clues to the nature of this enigmatic 
supergiant. 
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There are many mysteries regarding supernova physics including why a core collapse 
event produces an explosion releasing ~10 51 erg of kinetic energy and why some pulsars have 
space velocities as high as 1000 km s -1 (Strom, Johnston, Verbunt, & Aschenbach 1995). 
Both of these effects may be explained by anisotropics (Burrows 2000; Khokhlov et al. 1999; 
MacFadyen & Woosley 1999). Indeed, many core collapse supernova remnants show strong 
polarization (4% for SN1997X; Wang, Howell, Hoflich, & Wheeler 2001), indicating a high 
degree of asphericity. 

VY CMa exhibits asymmetries on scales from 10-10,000 R* (O'.'l - 100") and, therefore, 
is a good candidate for investigating the origin and nature of asymmetries in supergiants. 
Early on, Herbig (1972) showed that at optical wavelengths the nebula surrounding VY CMa 
appeared irregularly shaped with apparent structures within the bright central core of the 
nebulosity. At large scales (~ 3 pc), VY CMa's outflow is influencing the dynamics of the 
nearby H II region S310 (Lada & Reid 1978). Recent visible and near-infrared (near-IR) 
images of VY CMa from the Hubble Space Telescope (HST) display a complex distribution 
of knots and filaments embedded in a one-sided reflection nebula extending ~ 4" from the 
star (Smith et al. 2001; Kastner & Weintraub 1998) toward the southwest. On small scales, 
using near-IR aperture masking techniques, Monnier et al. (1999, 2004) have shown that the 
inner dust shell at ~ 15 R* is highly one-sided and consistent with the HST results. Even the 
molecular emission as traced by SiO is much brighter and more extended to the southwest 
(Shinnaga et al. 2003). Additionally, the spatial and redshift distribution of H 2 maser 
emission is not consistent with a simple spherical outflow, but rather appears elongated in 
the same direction as the near-IR emission (Marvel 1996; Richards, Yates, & Cohen 1998). 
However, the point-spread-function (PSF) subtracted 8.4 and 9.8 fim images of VY CMa by 
Smith et al. (2001) give the appearance of a bipolar structure within 2" east and west of the 
central source in apparent disagreement with the other observations. 

Compared to the mass-loss envelope's structure, little is known about the inner layers 
(< 10 R*) of VY CMa's atmosphere where the dust condenses and mass is driven out due 
to radiation pressure on the dust and subsequent drag on gas molecules. Since the star 
is concealed within an optically thick mass-loss envelope, only the longest wavelengths can 
possibly probe the inner atmosphere of VY CMa. Knapp, Bowers, Young, & Phillips (1995) 
found that the 8.4 GHz emission measured by the VLA is ~ 6 times that expected from 
photospheric emission of the supergiant and suggested that the excess may be due to an 
extended, 10,000 K ionized chromosphere as had been hypothesized for the nearby (D ~ 
130 pc) supergiant a Ori (Newell & Hjellming 1982; Hartmann & Avrett 1984). However, 
recent observations of a Ori revealed its atmosphere is dominated by cool (few x 1000 K) 
gas (Lim et al. 1998). Further, Lim et al. (1998), using the VLA in the A-configuration, 
found that a Ori was resolved at 43 GHz with a radius of ~ 44 mas and partially resolved at 
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8.4, 15, and 22 GHz. The radius at 43 GHz corresponds to a size of nearly twice the stellar 
radius of 22.5 mas (3 AU; Dyck, Benson, van Belle, & Ridgway 1996). Harper, Brown, & 
Lim (2001), using the flux density and visibility measurements from various authors and a 
modified version of the Mira radio photosphere model by (Reid & Menten 1997; hereafter 
RM97), created a semi-empirical model of the temperature and density distributions out to 
~ 14 R* that fit the radio measurements of a Ori . 

We have obtained VLA radio observations of VY CMa to test possible emission scenarios 
and probe for asymmetries at radio frequencies. We find that a radio photosphere model 
similar to RM97 and Harper, Brown, & Lim (2001) fit the observed data. We use this model 
to interpret the physical characteristics of the inner (1-10 R*) atmosphere of VY CMa. We 
also present new mid-IR observations of VY CMa at 11.7 and 17.9 fim and show that the 
thermal mid-IR emission from the envelope is consistent with that seen in optical and near-IR 
observations. 

This paper is organized as follows: in § 2 we describe the observations of VY CMa we 
have obtained with the VLA (§ 2.1) and the Keck Observatory (§ 2.2). We then discuss, in 
§ 3, the origin of the radio emission by first identifying emission models that do not fit the 
data well and then in § 3.1 present an atmospheric model containing a radio photosphere 
that reproduces the observed flux densities. A description of the analysis and discussion of 
our mid-IR data is contained in § 4. In § 5 we present our conclusions and suggest possible 
future observations. 



2. Observations 

2.1. VLA Radio Observations 

We observed VY CMa with the NRAO 1 Very Large Array (VLA) on 2001 Oct. 2 at 
43 GHz (Q-band) in the hybrid DnC configuration and on 2002 Sept. 21 and 23 at 22 GHz 
(K-band) and 15 GHz (U-band) in the CnB configuration; see Table 1 for details. We used 
the standard continuum observation technique, employing four 50 MHz bands (right- and 
left-circular polarization at two frequencies). The central frequencies observed were 43315 
and 43365 MHz, 22435 and 22485 MHz, and 14915 and 14965 MHz for Q-band, K-band, and 
U-band, respectively. Pointing was verified about every hour. For the 22 GHz and 15 GHz 
data, we employed fast-switching to remove the tropospheric phase fluctuations. The flux 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 



-4- 



density scale was set using observations of 3C 48 for the 43 GHz observations, and 3C 286 
for the 22 GHz and 15 GHz observations. The secondary calibrator, 0650 - 166, was used to 
calibrate the interferometric amplitudes and phases of VY CMa at all frequencies. Table 2 
lists the assumed calibration source flux densities. 

Calibration and imaging of the visibilities followed the procedures recommended for 
VLA data with the Astronomical Image Processing System (AIPS). Specifically, we imaged 
VY CMa at each wavelength with maps of 256 by 256 pixels and pixel scales of 0'.'14, O'/l, 
and O'/l at 43, 22, and 15 GHz, respectively. Due to the low declination of the source, 
we used hybrid VLA configurations with a a slightly longer northern arm which produced 
synthesized beams that are slightly elongated in the east-west direction: l'/7 x l'/5 with 
PA = -83° at 43 GHz; l'/O x0'/78 with PA = -70° at 22 GHz; and l'/5 x l'.'O with PA = 
80° at 15 GHz. At 43 and 15 GHz, VY CMa appears marginally larger than the beam, 
with deconvolved geometric-mean sizes of I'M and l'/7, respectively; at 22 GHz the source 
is unresolved. However, because the position angles for the elliptical Gaussian fits to the 
sources do not match those of the beam, we consider VY CMa to be unresolved at all three 
observed frequencies. 

We list the peak flux density measurements in Table 3, with formal errors that combine 
a 20% systematic error and the RMS noise error estimated from the background. We note 
that the integrated flux densities, from a Gaussian fit to the source, tend to be greater than 
the peak flux densities by as much as 20%. This is likely the result of poor phase calibration 
due to the > 10° separation between VY CMa and the secondary calibrator. We consider 
the peak flux density to be a better estimate of the flux since the source is unresolved at 
each wavelength. 

2.2. Keck Mid-IR Observations 

We observed VY CMa with the Long Wavelength Spectrograph (LWS; Jones & Puetter 
1993), a facility instrument at the W.M. Keck Observatory, on 2000 Feb. 5 at 11.7 fim (AA 
= 1.0 /im) and 17.9 fim (AA = 2.0 /im). LWS uses a 128 x 128 pixel Boeing Si:As detector 
and has a plate scale of 0'/08 pixel -1 , resulting in a 10'/2 x 10V2 field of view. We employed 
the "chop-nod" mode of observing, with a chop throw of 10" north. The IRTF IR standard 
a Aur (N = -1.94; Q = -1.93; Tokunaga 1984) was used to calibrate the flux of VY CMa and 
to measure the PSF. The full-width-at-half-maximum of a Aur was measured to be 0'.'58 at 
11.7 /im and 0748 at 17.9 /im. 

We measure fluxes of 8600±1700 Jy at 11.7 /im and 9300±1900 Jy at 17.9 ^m for VY 



- 5 - 



CMa. The 20% errors are derived from the variability of the calibrator's flux throughout the 
night. Our measured flux density at 11.7 /im is consistent with the color corrected (Beichman 
et al. f 988) IRAS 12 /im flux of 8729 Jy and with the mean DIRBE 12 /im flux of 9295.5 Jy. 

3. Identifying the Origin of the Radio Emission 

As we have not resolved VY CMa at radio wavelengths, we cannot comment on any 
asymmetry in the inner atmosphere, except to say that if one exists, it must be on scales 
smaller than ~ 1" = 1500 AU. The question we now ask is: what is the source of the radio 
emission from VY CMa? 

In Table 3, we have compiled the radio flux densities of VY CMa in the frequency range 
8.4 - 43 GHz from this work and the literature. From VLA monitoring of VY CMa at 8.4 
GHz from September 1990 to May 1991, Menten & Reid (2004) find that VY CMa may 
be variable by about a factor of 2 on an irregular timescale, but they note that within the 
errors, this variability is tentative. Knapp, Bowers, Young, & Phillips (1995) measured a flux 
density of 0.26 mJy for VY CMa in January 1991, consistent with the mean value of Menten 
& Reid (2004). Our measurement of VY CMa's flux density at 22 GHz of 0.66 ± 0.20 mJy is 
about 1/2 the 1.3 ± 0.3 mJy measured by Menten & Reid (2004) in the early 1990's at this 
frequency. This suggests that VY CMa has some significant variability at radio frequencies 
in addition to the irregular variability seen at optical wavelengths (Wallerstein 1978) and 
the semi-regular variability seen in maser emission (Harvey, Bechis, Wilson, & Ball 1974; 
Cox & Parker 1979; Gomez Balboa & Lepine 1986; Martinez, Bujarrabal, & Alcolea 1988). 
Further observations will be necessary to confirm any variability. 

Fitting a mean power law to the flux densities listed in Table 3 gives a spectral index 
of a ~ 1.4 ± 0.2 (Sj, oc u a ), but given the possible variability, this value is uncertain since 
the phase at each observation is unknown. If we break up the data into the '1990's set' and 
the '2000's set', where the former set is comprised of the Menten & Reid (2004) and Knapp, 
Bowers, Young, & Phillips (1995) data and the latter of the data presented herein, we find 
the spectral index of each set is 1.7 ± 0.3, which is similar to that found for the radio spectra 
of Mira variables (a = 1.86; RM97). This spectral index is also close that expected from 
the stellar photosphere (a — 2), however, the flux density levels we measure are about twice 
that expected from a star with photospheric temperature of T* = 2800 K and radius R* = 
2.25 x 10 14 cm (15 AU; Danchi et al. 1994) at the distance of VY CMa (see Figure 1, solid 
line) . 

The sub-millimeter spectrum of VY CMa was modeled by Knapp, Sandell, & Robson 
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1993; hereafter KSR93) using small circumstellar dust grains (a g = 2000 A; p g = 3.5 g cm -3 ) 
with an emissivity index f3 = 0.9. They derive a dust mass-loss rate of 2 x 10~ 6 M Q yr _1 
within a region of radius 10 18 cm. Assuming all the emission above the expected photospheric 
emission at radio frequencies is due to dust with the KSR93 parameters, then we calculate 
the total dust mass using 

^ _ 2(5 , „ - S phot )a g p g D 2 c 2 
3Q u kf ) T env h' 2 

where S v is the total measured flux density from the source, S p h ot is the expected photospheric 
flux density, D is the distance to VY CMa, Q u is the grain emissivity, and T env is the 
temperature in the envelope. We use the KSR93 function for the emissivity: Q u = Qo(v/voy 
where Qq = 5.65 x 10~ 4 at is = 274.6 GHz. The color corrected IRAS fluxes (Beichman 
et al. 1988) at 60 and 100 pm are 1119 Jy and 304.5 Jy, respectively, and the ratio of the 
IRAS fluxes gives an envelope temperature of 235 K. Thus, our 43.3 GHz flux density of 2.8 
mJy, after subtraction of 1.2 mJy from the expected photospheric emission, implies a dust 
mass of 1.3 x 10~ 3 M within our unresolved source of radius ~ Of 5. However, using the 
dust mass-loss rate and outflow velocity from the KSR93 model, the total mass expected to 
reside within a radius of 0V5 is only 1.8 x 10~ 4 M , a factor of ~ 7 less than that calculated 
from the 43.3 GHz data. That is, there would need to be more dust closer to the star than 
is predicted by the KSR93 model. It may be possible to fit the observed radio flux densities 
with a second population of dust grains closer to the star. 



3.1. Radio Photosphere Model 

We now turn to an alternative source of radio emission: a radio photosphere. In their 
attempt to model radio spectra of Mira stars, RM97 find that at long wavelengths their 
observations of these long period variables are well fit by radio photospheres with radii ~ 2 
R*. The opacity in the extended radio photosphere model is predominately provided by free- 
free absorption from the interaction between neutral atomic and molecular hydrogen and free 
electrons supplied mostly by Na, Al, K, and Ca (elements with low ionization potentials). 
We have applied this opacity model to VY CMa. 

Given a temperature and density profile for the stellar atmosphere as a function of radius 
outward from the surface, at each layer in the star's atmosphere we assume thermodynamic 
equilibrium and solve the Saha equation simultaneously for Na, Al, K, and Ca to determine 
the total electron density and pressure (P e ). Following RM97, we use solar abundances, 
ionization potentials, and partition functions from Gray (1992). We use dissociation coef- 
ficients given by Tsuji (1964) to calculate the densities of atomic and molecular hydrogen 
(p(H~) and p(H^)). To evaluate the absorption coefficients for H~ and H^~ free-free emission 
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(k u {R~) and /^(Hj)) for temperatures < 3000 K, we use the third order polynomials in tem- 
perature given in RM97. For temperatures higher than 3000 K, we assume the atmosphere 
is composed entirely of atomic hydrogen, and solve for the absorption coefficient, k„(H _ ), 
by numerically evaluating the equations given in Dalgarno & Lane (1966). The absorption 
coefficients are proportional to the inverse square of the frequency and have units of cm 4 
dyn -1 . We then calculate the optical depth increment 



where dl is the step size through the atmosphere. The brightness temperature, T b , at an 
optical depth r in the supergiant's atmosphere is computed from 



which is the formal solution to the radiative transfer equation in the Rayleigh- Jeans limit. 
We numerically integrate along the line-of-sight toward the center of the star and then repeat 
the integrations moving perpendicular to the line-of-sight with step size Sp to evaluate the 
center-to-limb brightness temperature profile, T b {p). The total flux density, S u , is calculated 
assuming azimuthal symmetry using 



p=0 

We use a model's total flux density values at all four observed radio frequencies compared 
to those observed to evaluate the validity of a model and its respective inputs result. 

As a simple approximation, we assume an isothermal atmosphere and compute the 
density structure for a range of input temperatures. The density profile for an isothermal 
atmosphere is that required by hydrostatic equilibrium out to the critical radius {r cr u = 
GeM*/2a 2 , where a = (k b T '/ ' jj) 1 ^ 2 is the isothermal sound speed, G is the gravitational 
constant, eM* is the effective stellar mass, and ji is the mean molecular weight of the gas) even 
if there is a radiation pressure driven wind (Lamers & Cassinelli 1999). We have used the term 
e = (1 — L^x/^cGM*), where x is the opacity, to account for the alteration of the effective 
stellar mass due to radiation pressure. We take the temperature- and density-dependent 
Planck mean opacity values from Helling, Winters, & Sedlmayr (2000) as upper limits to the 
opacity and consider this to be more accurate for the optically thin stellar atmosphere than 
the Rosseland mean opacities (Hoefner, Jorgensen, Loidl, & Aringer 1998). At temperatures 
higher than ~ 3000 K the Rosseland mean opacities are essentially equivalent to the Planck 
mean opacities (Helling, Winters, & Sedlmayr 2000). The density as a function of radius, r, 
for an atmosphere in hydrostatic equilibrium can be computed from 




(2) 




(3) 




(4) 




(5) 
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where p(R*) is the density at R*, and H* = k b T*/fig* is the scale height with = GeM^/R 2 . 

We use a stellar radius of R* = 2.25 x 10 14 cm, an initial density of = 5 x 10 14 cm -3 , 
and a stellar mass of M* = 15 M . The initial density has been estimated by assuming a 
2800 K optical photosphere in hydrostatic equilibrium with a Rosseland mean opacity of x 
= 10" 4 cm 2 g _1 at an optical depth of r = 2/3. We report the results of three isothermal 
model atmospheres: T = 1600 K, 3000 K, and 4500 K. Figure 2 shows the density profile for 
each model while Figure 3 shows the output brightness temperature as a function of impact 
parameter for each model at 22 GHz. 

For the T = 1600 K model atmosphere, the density drops so quickly that the optical 
depth toward impact parameters > R* is very small and, therefore, the radius of the radio 
star is not much larger than R*. Since is a sum over the radio radius, this results in flux 
density values too low to be consistent with those measured for VY CMa. 

We find that the two warmer atmospheres produce radio photospheres whose model 
flux densities bracket those observed for VY CMa. The model output flux density values 
from the radio photosphere for both these models are listed in Table 3 and are represented 
as the dashed and dot-dashed lines in Figure 1. The spectral index for both models is 1.8. 
Because of the optically thick envelope enshrouding VY CMa, the supergiant lacks any UV 
lines which would hint at the existence and temperature of a stellar chromosphere. Recent 
modeling of high resolution UV spectra of the supergiant a Ori's atmosphere produced a 
temperature fit to the Ha lines not in excess of 5500 K (Lobel 2003). By analogy, such a 
distribution of warm plasma may exist around VY CMa. 

Using these models, we find that the apparent radius at radio frequencies, estimated as 
the half-width at half-brightness of the T b (p) profile, is 1.5-2 R* and varies inversely with 
frequency; see Figure 4. The dependence of apparent radius on frequency is consistent with 
the pattern observed by Lim et al. (1998) for the resolved and partially resolved radii of a 
Ori. The apparent radius of VY CMa's model radio photosphere is comparable to the 2 R* 
radio photosphere proposed for Mira stars (RM97) and is consistent with the radio emitting 
region fit by semi-empirical models in Harper, Brown, & Lim (2001) for a Ori (see their 
Figure 5). 

The mass-loss rate for an isothermal atmosphere can be calculated from 

(corrected from Lamers & Cassinelli 1999), where v esc (R*) = (2 G eM*//?*) 1 / 2 is the escape 
velocity at the stellar photosphere. For the T = 3000 K model, then a = 5.0 km s -1 and 
v esc (R*) = 42.1 km s" 1 , which gives M tot = 1.9 x 10 -12 M Q yr -1 , about 8 orders of magnitude 



- 9- 



lower than what has been measured for VY CMa. For the T = 4500 K model, however, a = 
6.1 km s _1 and v esc (R*) = 33.9 km s _1 , which gives M tot = 2.8 x 10~ 4 M yr" 1 , comparable 
to the value measured by molecular line emission of 2-4 x 10~ 4 M Q yr _1 (see Danchi et al. 
(1994) for a summary). 

We have also applied the radio photosphere model to an adiabatic atmosphere. How- 
ever, with an adiabatic atmosphere in hydrostatic equilibrium, we find that the optical depth 
due to the H~ and H^~ opacity never reaches > 1. At high temperatures, where the Planck 
mean opacity is large (x = 0.32 for T* = 6000 K), for the atmosphere to be in hydrostatic 
equilibrium the density as a function of radius must drop faster compared to lower temper- 
atures. Once the density falls below a threshold density (n# ~ 10 12 cm -3 ), there simply are 
not enough free electrons for the H~ and H2 opacity to be effective. 



4. Mid-IR Emission 

Our new resolved mid-IR images of VY CMa with the LWS at the W. M. Keck Ob- 
servatory are shown in Figures 5a & c. We find that VY CMa is asymmetrically extended 
to the southwest at both 11.7 and 17.9 /im with the maximum extension at position angle 
277° as measured East from North. Figure 5b & d show the radial profiles at maximum 
extension (277°), minimum extension (47°), and at the normal to the maximum extension 
(137° & 317°) compared to the azimuthally averaged radial profile of the standard star a 
Aur (solid line with error bars). We find that at 11.7 (17.9) //m, at a distance of Of! 8, the 
flux is 3.4 (2.0) times greater toward position angle 227° compared to along position angle 
47°. At a distance of l'.'O, the ratio of fluxes is 3.4 (3.5) at 11.7 (17.9) /xm. The direction 
of the extended emission is consistent with that seen in the optical and near-IR. We find no 
evidence for the bi-polar east-west structure surrounding the central star at a distance of ~ 
2" as seen in the PSF-subtracted images of Smith et al. (2001). The average polarization of 
the SiO masers in the v — 1, J = 5-4 transition, 251°, is essentially aligned with the max- 
imum extent of the mid-IR emission (Shinnaga, Moran, Young & Ho 2003). This implies 
that the asymmetrical envelope is not just caused by a scattering effect, but that the mass 
itself is anisotropically distributed, perhaps by a bipolar outflow. 

The color temperature of the dust grains is estimated to be ~ 245 K from the ratio 
F nu (11.7 /im)/F nu (17.9 /im) = 0.92. We estimate the dust mass M dust using 

F U D 2 

M dust = r> Irp \ ( 7 ) 

where F v is the observed flux at 17.9 jum and \v is the grain opacity. With T d = 245K, and 
X„(17.9 /im) ~ 600 cm g" 1 (Sopka et al. 1985), we find M dust = 5.8 x 10~ 4 M . This is a lower 
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limit to the dust mass since we have assumed the circumstellar envelope is optically thin in 
this calculation. The dust mass-loss rate (M dust ) can be estimated using the approximate 
size of the dust emitting region (r dust ~ 2") and the outflow velocity (v out = 39 km s _1 ; 
Sopka et al. 1985) and is found to be 1.6 x 10~ 6 M yr^ 1 , which is consistent with the value 
obtained from 400 /zm data by Sopka et al. (1985) of 7.8 x 10~ 7 M Q yr _1 and the upper 
limit of 3.5 x 10~ 6 M Q yr _1 published by KSR93 and calculated from 1.1 mm data. 

5. Conclusions 

We report the detection of emission from the supergiant VY CMa at radio frequencies 
with the VLA. We have constructed isothermal outer atmospheres of 3000-4500 K containing 
radio photospheres where interactions between free electrons and neutral and molecular 
hydrogen provide the dominant opacity. These isothermal atmospheric models produce a 
radio photosphere at 1.5-2 R* and flux density values and a spectral index comparable to 
those measured for VY CMa. An extrapolation of the model can account for the observed 
total mass-loss rate observed for VY CMa. We also present resolved mid-IR imaging of VY 
CMa. These data suggest that warm dust is extended in the same direction as the near-IR 
reflection nebula around VY CMa and that the mass-loss rate as derived from the mid-IR 
is consistent with previous measurements. 

In their VLA observations of a Ori, Lim et al. (1998) resolved the supergiant's radio 
surface at 43 GHz and found that the star was asymmetrically shaped at this frequency. The 
authors proposed that such an asymmetry may arise from giant convection cells elevating 
cool material into the atmosphere. Our VLA radio observations at 15, 22, and 43 GHz 
produce unresolved images of VY CMa and thus we cannot comment on the asymmetry 
in the supergiant's inner atmosphere. To resolve a radio photosphere at 1.5-2 R* would 
require a resolution of ~ 0'/04 which is just beyond the capability of the VLA at 43 GHz. 
However, it will be possible to resolve VY CMa at higher frequencies with the Atacama Large 
Millimeter Array (ALMA). We expect a flux density at 345 GHz of ~ 120-240 mJy with an 
apparent radius of ~ 0'.'025, which would be easily resolved with ALMA's largest D = 10 km 
configuration. A resolved detection would confirm the existence of a radio photosphere in VY 
CMa and may provide clues to the symmetry of such an extended photosphere. Continued 
monitoring of VY CMa at radio frequencies will help identify any variability. Further, 
monitoring of VY CMa in the IR as proposed by Monnier et al. (2004) and Smith et al. 
(2001) to follow evolution of structures and the general shape of the mass-loss envelope will 
provide insight into the orientation and dynamics of the central star and may help constrain 
the dust mass-loss parameters. 
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Fig. 1. — Measured flux densities of VY CMa from 8.4 to 43 GHz from Table 3. The 
solid line represents a blackbody curve for a 2800K star with radius 2.25 x 10 13 cm. The 
dashed (dotted) line connects the output flux densities from the 3000 K (4500 K) isothermal 
atmospheric model with a radio photosphere (see § 3.1 for details). 
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Table 1. VLA Observations of VY CMa 



Date 


Freq. 


Config. 


On-Source Time 


Beam Size 


PA 




(GHz) 




(hr) 




(deg) 


2001 Oct. 2 


43 


DnC 


3.2 


1"7 x r?5 


-83 


2002 Sept. 21 & 23 


22 


CnB 


3.3 


l'.'O x 0"8 


-70 


2002 Sept. 21 & 23 


15 


CnB 


3.3 


1'.'5 x r?2 


80 
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Table 2. Calibrators 



Source 


R.A. 


Dec. 




^22GHz 






(J2000) 


(J2000) 


(Jy) 


(Jy) 


(Jy) 


3C48 


01 37 41.2999 


+33 09 35.133 






0.53 


3C286 


13 31 08.2880 


+30 30 32.959 


3.45 


2.52 




0650-166 


06 50 24.5819 


-16 37 39.725 


2.55 


2.38 


1.18 



Table 3. Observed and Model Flux Densities for VY CMa 



Frequency 
(GHz) 


Measured Y v 
(mJy) 


F„(T = 3000 K Model) 
(mJy) 


F„(T = 4500 K Model) 
(mJy) 


8.4 


0.26 ± 0.03 a ; 0.24 ± 0.02 b 


0.10 


0.24 


15 


0.41 ± 0.092 c 


0.29 


0.70 


22 


1.3 ± 0.3 b ; 0.69 ± 0.15 c 


0.61 


1.4 


43 


2.5 ± 0.51 c 


2.2 


4.9 



a Knapp, Bowers, Young, & Phillips (1995); data from January 1991 
b Menten & Reid (2004); data from September 1990 to May 1991 
This work 
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Fig. 2. — The density profiles for isothermal stellar atmospheres with R* = 2.25 x 10 14 cm, 
p, = 5x 10 14 cm" 3 , and M* = 15 M . Dot-dashed line: T = 1600 K; dashed line: T = 
3000 K; dotted line: T = 4500 K. 

Fig. 3. — Calculated brightness temperature from H~ and free- free interactions at a 
frequency of 22 GHz for the isothermal stellar atmosphere models in Figure 2. 

Fig. 4. — Apparent stellar radius at radio wavelengths as a function of frequency for both 
3000 K (dashed line) and 4500 K (dotted line) isothermal radio photosphere models described 
in the paper and depicted in Figures 2 & 3. 

Fig. 5. — a) The 11.7 /zm image of VY CMa. North is up and East to the left, b) Cuts through 
the 11.7 /iin image at angles representing the semi-major axis of the extended source (227° 
and 47°) and the semi-minor axis (137° and 317°) compared to the azimuthally averaged 
PSF of the standard star a Aur (solid line with error bars). The angles are measured East 
from North, c) the 17.9 /zm image of VY CMa. d) same as b) for the 17.9 /iin data. 
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